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Motivated by developments in quantum information science, much recent effort has been directed 
toward coupling individual quantum emitters to optical microcavities. Such systems can be used 
0^ ' to produce single photons on demand, enable nonlinear optical switching at a single photon level, 

' and implement functional nodes of a quantum network, where the emitters serve as processing 

nodes and photons are used for long-distance quantum communication. For many of these practical 
applications, it is important to develop techniques that allow one to generate outgoing single photons 

' of desired frequency and bandwidth, enabling hybrid networks connecting different types of emitters 

and long-distance transmission over telecommunications wavelengths. Here, we propose a novel 
approach that makes use of a nonlinear optical resonator, in which the single photon originating 
from the atom-like emitter is directly converted into a photon with desired frequency and bandwidth 

\ using the intracavity nonlinearity. As specific examples, we discuss a high-finesse, TE-TM double- 

mode photonic crystal cavity design that allows for direct generation of single photons at telecom 
wavelengths starting from an InAs/GaAs quantum dot with a 950 nm transition wavelength, and a 
-h | scheme for direct optical coupling of such a quantum dot with a diamond nitrogen-vacancy center 

at 637 nm. 

Nonlinear optical frequency conversion is widely used in fields as diverse as ultrahigh-resolution imaging M and 
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telecommunications, as it allows for the generation of light in parts of the spectrum for which there are no convenient 
sources. Common implementations include optical parametric oscillators to make tunable femtosecond lasers in the 
infrared, and conversion of the 1064 nm Nd:YAG laser to make green laser sources via second-harmonic generation. 
Recently, nonlinear photonic crystal cavities @> S S @] have emerged as promising systems in which similar 
nonlinear functionalities can be achieved at micron scales, which would enable the miniaturization of optical devices 
onto integrated platforms. While the majority of such work focuses on conversion of classical fields, these systems are 
now being applied to quantum optics and quantum information science @, H, H, E3| • 
^ In recent years, there has also been a concerted research effort to develop on-demand single- pho ton sources using 

single quantum emitters strongly coupled to resonant optical microcavities (cavity QED) [TJ 

Hi El- The str one 

| coupling of the emitter to a resonant cavity results in preferential emission into the cavity mode of a single photon 
with frequency near the atomic resonance. Connecting pairs of such systems would form the basis for distributed 
quantum networks, where the emitters serve as processors and photons carry information between the nodes [l4T |. In 
practice, however, the photon emission occurs at wavelengths determined by the atomic resonance frequency. This is 
impractical, as it does not exploit the low-loss telecom frequency band for long-distance transmission and requires all 
• i-h , emitters in a quantum network to be identical. 

Here, we describe a novel approach to generate single photons with controllable wavelength and bandwidth. Our 
approach makes use of an integrated nonlinear optical cavity in which optical emission is directly frequency-shifted 
into the desired domain using intracavity nonlinear optical processes. This cavity-based generation technique is quite 
robust in that the maximum efficiency does not depend on an explicit phase-matching condition [l5j , as would occur in 
an extended nonlinear crystal or fiber, but rather only on the ratio of the cavity quality factor to mode volume (Q/V). 
As an example, we demonstrate a novel double mode TE-TM cavity design in a GaAs photonic crystal that is well- 
suited for the conversion of photons from quantum dots to the telecom band. We also present a similar GaP-based 
desi gn f or direct coupling between a nitrogen-vacancy center in diamond [l6l UtI [Lsl. Il9j and an InAs/GaAs quantum 
dot [20l . l2ll . |22| , which could enable practical realization of a heterogeneous quantum network. In addition to effective 
wavelength control of single photons @, H, S US] , our approach enables the manipulation of their bandwidth, which 
is important for fast communication. 
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The concept of single-photon spectral control 

We first discuss the general protocol for generating single photons on demand at arbitrary frequencies using a 
nonlinear double-mode cavity, and introduce a simple theoretical model to derive the efficiency of the process. The 
system of interest is illustrated schematically in Fig. [TJ As in standard cavity-based single-photon generation proto- 
cols [l3l [l4| , a single three-level atom (or any other quantum dipole emitter) is resonantly coupled to one mode (here 
denoted a) of an optical cavity. The emitter is initialized in metastable state |s), and an external laser field with 
controllable Rabi frequency fi(i) couples |s) to excited state |e). The transition between |e) and ground state \g) 
is resonantly coupled to cavity mode a (frequency uj a ), with a single-photon Rabi frequency g\. The relevant decay 
mechanisms (illustrated with gray arrows in the figure) are a leakage rate n a for photons to leave cavity mode a, and 
a rate 7 that |e) spontaneously emits into free space rather than into the cavity. Conventionally, in absence of an 
optical nonlinearity, the control field Q(t) creates a single atomic excitation at some desired time in the system, which 
via the coupling g\ is converted into a single, resonant cavity photon. This photon eventually leaks out of the cavity 
and constitutes an outgoing, resonant single photon generated on demand whose spatial wave-packet can be shaped 
by properly choosing Q(t) fl4j ]. 

In our system, the cavity is also assumed to possess a second mode c with frequency u) c , and our goal is to induce 
the single photon to exit at this frequency rather than uj a . This can be achieved, provided that the cavity medium 
itself possesses a second-order (x^) nonlinear susceptibility, by applying a classical pump field to the system at the 
difference frequency tub = t-u a — uj c . The induced coherent coupling rate between modes a and c is denoted gi- The 
field b need not correspond to a cavity mode. Mode c has a photon leakage rate, which we separate into an "inherent" 
rate, K Ct i n , and a "desirable" (extrinsic) rate, K c ^ ex . K c .in characterizes the natural leakage into radiation modes and 
also absorption losses, and can be expressed in terms of the (unloaded) cavity quality factor as K c .in = ^c/^Qc- K c ^ ex 
characterizes the out-coupling rate into any external waveguide used for photon extraction. The total leakage of mode 
c is then k c = K c ^ n + n c ^ ex . 

More quantitatively, the effective Hamiltonian for the system (in a rotating frame) is given by 

Hi = H c + Hi oss , 

H c = h,gi(a eg a a + a ge al) + h,n(t)(a es +a se ) + hg 2 (ala c + a a al), 
Hi oas = -ycr ee - — a> a a a ■ — ■ — a l c a c , (1) 

where H c describes the coherent part of the system evolution (for simplicity we take gi 2, to be real), and Hi oss is 
a non-Hermitian term characterizing the losses. <xy = are atomic operators, while ai is the photon annihilation 



a) b) |e,O a A) 




FIG. 1: Schematic of single-photon frequency conversion a) A single three-level emitter is coupled to a double-mode 
cavity that possesses a x' 2 ' nonlinearity. After excitation, the emitter emits a photon into the cavity at frequency w a . When 
the cavity is irradiated by the pump beam at cot, the photon is converted to a second cavity mode at frequency uj c . b) Level 
diagram: coherent coupling strengths are indicated with blue arrows, while gray arrows denote undesirable loss mechanisms. 
The emitter is controllably pumped from initial state \s) via an external laser field Q.(t) to excited state |e). The excited state 
|e) can reversibly emit a single photon into cavity mode a (while bringing the atom into state |<?)) at a rate g\, and can also 
decay into free space at rate 7. Mode a has an inherent decay rate given by K a . The nonlinearity allows the photon in mode 
a to be converted to one in mode c at a rate 52 when the cavity is pumped by a laser of frequency ujb — u a — ui c . The leakage 
rate of the frequency-converted photon at ui c is split up into undesirable channels (K c ,in) and desirable out-coupling to a nearby 
waveguide (k c ^x)- 
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operator for mode i. The vacuum Rabi splitting gi can be written in the form gi = d • E a (r)/7i, where d is the dipole 
matrix element of the |<?)-|e) transition, and E a (r) is the electric field amplitude per photon at the emitter position 
r. The electric field per photon in mode i = a, c is determined by the normalization 

^ = J dr £oe (r)|E 4 (r)| 2 , (2) 

where e(r) is the dimensionless electric permittivity of the material. The nonlinearity parameter is given by [|[ 

.92 = -~ J dvxfXKi (E btj E c , k + E Ctj E btk ) . (3) 

The amplitudes E a ^ c appearing above are normalized by Eq. while Ef, is the classical pump amplitude. Importantly, 
one can compensate for a small nonlinear susceptibility x or held overlap (phase matching) simply by using larger 
pump amplitudes E}, to achieve a desired g2 strength. 

For a system initialized in |s), there can never be more than one excitation, and the system generally exists as 
a superposition of having the system in state \s) or |e) (with no photons) or having a photon in one of the modes 
a,c (and the emitter in |g)), 

\ip(t)) = c s {t)\s) + Ce(t)\e) + c a (t)\l a ) + c c {t)\l c ). (4) 

The system is initialized to c s (0) = 1 with all other Cj(0) = and the time evolution is given by Cj = —(i/Ti)(j\Hi\il>(t)). 
In this effective wave- function approach, provided that \s) is always driven, Y] ^ \cj\ —*0 as t^oo due to losses, which 
can be connected with population leakage out of one of the aforementioned decay channels. In the limit that f2(t) is 
small and varies slowly, all other Ci(t) adiabatically follow c s (t) (see the Methods section), and one finds 

c e (t)*-m(t)c s (t)-l (i+ — 4 f\, )c e . (5) 



2 V «a + 4g 2 /K c 

Physically, we can identify 7totai = 7 + ~ +lg*/ K as the cavity-enhanced total decay rate of |e), where the first (sec- 
ond) term corresponds to direct radiative emission (emission into mode a). Similarly, the denominator n a + \g\j ' k c 
corresponds to the new total decay rate of mode a in the presence of an optical nonlinearity, as it yields a new channel 
for photons to effectively "decay" out of mode a into c at rate \g\j k c . It is clear that some optimal value of g^ 
exists for frequency conversion to occur. In particular, for no nonlinearity (52 = 0) this probability is non-existent. 
On the other hand, for gi — > 00, one finds 7totai = 7, which indicates that the leakage from mode a into c is so 
strong that no cavity-enhanced emission occurs. Note that the use of time-varying control and pump fields allows 
for arbitrary shaping of the outgoing single-photon wavepacket at frequency u> c , provided only that the photon band- 
width is smaller than k c (physically, the photon cannot leave faster than the rate determined by the cavity decay, 
see Methods). This feature is particularly useful in two respects. First, in practice k c can be much larger than 7, 
which enables extremely fast operation times. Second, pulse shaping is useful for constructing quantum networks, as 
it allows one to impedance-match the outgoing photon to other nodes of the network. 

Based on the above arguments, the probability that a single photon of frequency lo c is produced and extracted into 
the desired out-coupling waveguide is given by 

p Ci n <j) K c ex 

~ 1 + (j) + C in 1 + (f> K C ' 

where tf> — Ag^ / {n a K c ) characterizes the branching ratio in mode a of nonlinearity-induced leakage to inherent losses, 
and Ci n = 4<? 2 /7K a is the inherent cavity cooperativity parameter for mode a in absence of nonlinearity. The first 
term on the right denotes the probability for |e) to decay into mode a, the second term the probability that a photon 
in mode a couples into mode c, and the third term the probability that a photon in mode c out-couples into the 
desired channel (see Methods for an exact calculation). <j) depends on the pump amplitude Eb, with the optimal value 
4> = \/l + Ci n yielding the maximum in F. For large Ci„3>l, the maximum probability is 

Fa(l--iLJ)^. (7) 



Considering an emitter placed near the field maximum of mode a, Cin^^S- n &y —, where Q a , V a are the mode quality 
factor and volume, respectively, and n is the index of refraction at frequency uo a . The ratio 7/70 is the spontaneous 
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FIG. 2: Cavity mode characteristics Frequency conversion platform based on a photonic crystal nanobeam cavity, integrated 
extraction waveguide, and off-chip coupling laser (uib) tuned to the difference frequency of the modes. The cavity is formed by 
introducing a local perturbation into a periodic ID line of air holes in the free-standing nanobeam. The desirable (« Ciea; ) and 
inherent (/t c ,in) loss channels from mode c are shown. The insets show the schematic cavity spectrum with photonic stopbands 
shown in grey, and the dominant field components of the TEo and TM2 (w a ) modes. The t/z-plane cross-sections of the 
modes (upper left) show the E y (E z ) component of mode c (a) at the center of the cavity, highlighting the mode overlap and 
polarizations. In the optimized structure (described in Methods), the TE mode at 1425 nm has Q = 1.2 x 10 7 and V n = 0.77, 
and the TM mode at 950 nm has Q — 7.3 x 10 4 and V„ = 1.45 (V n is the mode volume normalized by (A/n) 3 ). The inherent 
peak cooperativities for the modes are C; n = 2.4 x 10 and C in = 3.7 x 10 , which are well into the strong coupling regime, 
as given by C > 1. 



emission rate into non-cavity modes normalized by the spontaneous emission rate 70=710^ |d| 2 / (37reo?ic 3 ) of an emitter 
embedded in an isotropic medium of index n. This ratio is expected to be of order 1 — 10 for our devices of interest, 
and thus the efficiency essentially depends only on Q a /V a . 

Finally, while we have focused on the case of single-photon generation here, the reverse process can also be con- 
sidered, where a single incoming photon at frequency lo c is incident upon the system, converted into a photon in 
mode a, and coherently absorbed by an atom with the aid of an impedance-matched pulse Sl(i), causing its internal 
state to flip from \g) to \s). Generally, by time-reversal arguments [23], it can be shown that the probability F for 
single-photon storage is the same as that for generation. 



Realization in a nonlinear photonic crystal cavity 



In order to implement this frequency conversion scheme in a practical fashion, there are several constraints on the 
design of the cavity modes. For the nonlinear process to be efficient, mode a must have a high cooperativity (Q/V) to 
ensure strong coupling of the emitter (see Fig.[T]). For mode c, a high Q factor (small n c ) is important to maximize the 
nonlinear coupling parameter, </>, and hence reduce the pump power needed in order to reach the optimum nonlinear 
coupling strength, gi- The cavity should also be composed of a x^ nonlinear material that is transparent in the 
desired frequency range. Finally, in order for the modes to couple efficiently via the nonlinear susceptibility of the 
cavity, they must have a large spatial overlap and the appropriate vector orientation, as determined by the elements 
of the x^ 2 ) tensor of the cavity material (see Eq. ([3])). 

As a host platform for the nonlinear cavity, the III-V semiconductors are promising candidates because of their 
significant second-order nonlinear susceptibilities and mature nanofabrication technologies. However, the symmetry of 
the III-V group x^ tensor (x^l-T^O, * 7^.7 7^) requires that the dominant field components of the modes be orthogonal 
in order to maximize the nonlinear coupling. It further implies that if the classical field which drives the nonlinear 
polarization is incident from the normal direction (e.g., from an off-chip laser), one of the cavity modes must have a 
TM polarization. 
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We adopt a photonic crystal platform to realize a wavelength-scale nonlinear cavity that meets these requirements. 
Recently, 2D photonic crystal nanocavities have shown great promise for strongly coupling an optical mode to a 
quantum dot emitter [2l|, [22j]. In addition, they have been used as platforms for classical nonlinear optical generation 
and switching @)H3]- The challenge, however, is to design a nonlinear photonic crystal nanocavity which supports 
two orthogonal, high cooperativity modes with a large mode field overlap. 

To enable a monolithic cavity design which supports both TE and TM modes, we design a photonic crystal 
"nanobeam" cavity - a free-standing ridge waveguide patterned with a one-dimensional (ID) lattice of holes - for 
which we can control both TE and TM photonic bandstructures. Recently, there has been much interest in photonic 
crystal nanobeam cavities [HI, [H, [2?], HH, Hi| due to their exceptional cavity figures of merit (Q and V) , relative ease 
of design and fabrication, and potential as a platform to realize novel optomechanical effects (30l . l3l| . Our frequency 
conversion scheme can be realized in a similar structure, as shown in Fig. [51 We optimize two high cooperativity 
cavity modes by exploiting the different quasi-lD TE and TM photonic stopbands of the patterned nanobeam (shaded 
regions in the inset of Fig. [2]). A key design point is that the TE and TM bandstructures can be tuned somewhat 
independently by varying the cross-sectional aspect ratio of the ridge. For example, in a nanobeam with a square 
cross-section, the two stopbands overlap completely. As the width-to-depth ratio of the waveguide is increased, the 
effective index of the TE modes increases relative to the TM modes, shifting the TE stopband to longer wavelengths. 

Example implementations 

As a first example, we design a GaAs photonic crystal nanobeam cavity with modes at 950 nm and 1425 nm suitable 
for directly generating single photons at telecom wavelengths from InAs/GaAs quantum dots. To achieve such a large 
spectral separation, we couple the fundamental TE cavity mode to a higher-order TM 2 cavity mode (see inset Fig. [5]). 
Crucially, the photonic crystal lattice tapering [2?], [H, Hl[ is effective in enhancing the Q factors of both TE and TM 
modes. Details of the cavity parameters and optimization are provided in the Methods section. For this cavity, the 
coupling field (Pt,) must have a wavelength A;, = 2.85 /im in order to efficiently drive the difference- frequency process. 
GaAs is an attractive nonlinear cavity material because it has a reasonably large strength [32|], a high refractive 
index, and mature microfabrication techniques. 

As evident in Fig. [2 the overlap of the two modes changes sign near the edges of the ridge compared to the 
middle due to the different symmetries of the TEo and TM2 modes. However, the induced nonlinear polarization is 
dominated by the negatively signed anti-nodes near the middle of the ridge, and the imperfect overlap in the integral 
can be completely compensated for by a stronger pump beam. Thus, by selecting a higher order TM2 mode, we 
have gained a larger frequency conversion bandwidth at the expense of the somewhat higher pump power required to 
overcome the ensuing phase mismatch. Note that the fundamental TEq-TMo mode overlap is nearly ideal, and would 
be appropriate for applications requiring relatively small frequency shifts. 

We now calculate the probability to convert a single photon from 950 to 1425 nm in our system. The optimized 
cavity design (see Methods) simultaneously yields high quality factors and small mode volumes, which allows for 
extremely high cooperativities for each mode (C„ > 10 4 ). From Eq. Q we find that this enables an internal 
conversion probability of up to F — 0.99 when waveguide extraction efficiency is not taken into account. In practice, 
to efficiently out-couple the frequency-converted single photon into a waveguide, we require the ratio S = K c ^ ex / 'n c .m 
to be large (i.e. overcouplcd). The branching parameter <j) scales as p,/ (K c ,m(l + S)), and so to increase the extraction 
ratio <5, the pump power (p,) must also be increased to maintain the optimal <fi. Essentially, achieving good extraction 
efficiency requires one to intentionally increase the losses in mode c (via the out-coupling waveguide), which in turn 
requires more pump power to maintain the critical coupling. This relationship is made clear in Fig. [3] which plots 
the probability F as a function of pump power p, and extraction ratio 5. For a given 5, the power F& can be chosen 
to maximize the probability, reflecting the optimal value of g<i for frequency conversion. The probability rises rapidly 
with p,, reaching a maximum at relatively low powers (visible as the sharp ridge in the contours). Three fixed 
S contours are plotted in Fig. [3]Jb), demonstrating that efficient extraction of frequency-converted photons can be 
realized at modest pumping powers. For example, for 8 — 10, an extraction probability of 0.7 can be realized with 
a coupling laser power of 3 mW focused in a diffraction-limited focal spot. We note that the absorption of GaAs at 
2.85 fiia is negligible, and so there will be no pump-induced heating. For this particular cavity design, the outgoing 
converted photon can be shaped to have a bandwidth of up to k c ~100 MHz. 

By exploiting the scaling properties of Maxwell's equations, it is straightforward to design a similar cavity in GaP 
which supports modes at 637 nm and 950 nm. GaP is a nonlinear wide bandgap semiconductor which is transparent at 
637 nm and has recently shown promise for the microcavity enhancement of diamond NV emission [33l . |34| . Accounting 
for the exact refractive index dispersion and x*- 2 ^ strength of GaP, we calculate that the internal frequency conversion 
probability is 0.99, and the extraction probability is 0.7 for < 4 mW coupling power. The 637-950 nm span would be 
sufficient to couple any pair of the most relevant quantum emitters, namely NV centers in diamond; atoms such as 
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(a) (b) 

FIG. 3: Probability of single-photon frequency conversion from 950 nm to 1425 nm. The photon is coupled into a 
well-defined output channel at rate n c ^ x . Note that the internal probability of conversion in the absence of an over-coupled 
extraction channel is 0.99. (a) Probability as a function of the pump laser power, Pb, and the extraction ratio, 5 = K CtBX / K C) i n . 
For a given 5, there is an optimal operating power, Pb, as visible by the sharp contour ridge at small Pb. (b) Probability as 
a function of Pb for different values of <5. Because of the rapid rise in probability at low Pi,, the system does not need to be 
operated at the optimum to achieve high conversion probabilities. For example, for 8 — 10, a probability of F = 0.7 can be 
achieved with a pump power Pb = 3 m W (indicated by the arrow) . 

Cs or Rb; and InAs/GaAs quantum dots. Such a cavity could also be integral to creating a stable, room temperature 
single-photon source emitting in the telecom band based on frequency-converted NV center emission in diamond [35| . 
Given that it may be difficult to span the large spectrum from 637 nm to telecom wavelengths in a single monolithic 
design, a 637-950 nm cavity could be the first stage of a two-step frequency conversion process involving our first 
example as the second stage. More generally, cascading allows our design to be extended to cover virtually any 
frequency span. 

Outlook 

We have shown that high-fidelity, intra-cavity frequency conversion of single photons from a dipole-like emitter 
can be achieved using a two-mode nonlinear cavity pumped by a classical field. Our general framework is valid for 
conversion between arbitrary frequencies, and the efficiency depends only on the cavity parameter Q/V. As realistic 
implementations, we propose two different high-cooperativity, double-mode photonic crystal nanocavities to enable 
highly efficient coupling from 950 nm - 1425 nm and 637 nm - 950 nm, respectively. Single-photon conversion between 
these wavelengths would allow, in the first instance, coupling of InAs/GaAs quantum dot emission [2fj, l2l|, [22J into low- 
loss optical fibers in the telecom spectrum. The second example would facilitate a direct optical connection between 
two types of solid-state emitters currently of great interest: a nitrogen- vacancy center in diamond [IB, Ell Ell, Ell and 
an InAs/GaAs quantum dot. Integrated together, the two designs could allow for cascaded frequency conversion of 
NV center emission into the telecom band. Further design improvements should lead to larger frequency spans and 
also lower pump power requirements (e.g., by allowing LOh to correspond to a third cavity mode). Although we have 
emphasized large frequency shifts in this paper, a smaller shift could be readily achieved by coupling the TEo mode 
with the fundamental TMo mode, which has a larger Q factor than the TM2 mode studied here. The TE0-TM0 modes 
have a larger spatial overlap, reducing the coupling power required for high probability frequency conversion. 

Beyond the aforementioned applications, the techniques described here can potentially be extended to open up many 
intriguing opportunities. For example, the photon emission of a particular emitter could be shifted into wavelengths 
where high-efficiency detectors are available. It also allows coupling of atomic emitters such as Cs or Rb with solid- 
state emitters to create hybrid atom-photonic chips f36{ . In addition, a number of quantum entanglement schemes 
for atoms rely on joint photon emission and subsequent detection to probabilistically project the atomic system into 
an entangled state [37], [38|, l39( . Such schemes rely on the indistinguishability of photons emitted from each atom, and 
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implementing such techniques in nonlinear cavities could allow entanglement between different types of emitters. In 
addition, the protocol described here could be extended for generating narrow-bandwidth, entangled photon pairs with 
high efficiency and repetition rates, which are a valuable resource for applications such as quantum cryptography [Zol ]. 
Our scheme could also be applicable in active materials, where laser wavelengths could be converted from easily 
accessible regions like 1500 nm to the mid-infrared range. Finally, it would be intriguing to combine these ideas 
with cavities that exhibit opto-mechanical coupling (3lj | , which would potentially allow the photonic frequencies to be 
dynamically and rapidly tuned. 



Methods 



Derivation of nonlinear conversion efficiency 

The state amplitudes of the wave-function given in Eq. @ evolve under the interaction Hamiltonian Hj of Eq. (TT]) 
through the following equations, 

c s = -iQ(t)c e , 

c e = ~iQ(t)c s - igic a - ( r r/2)c e , 

Ca = -igiCe - ig2C c ~ (K a /2)c a , 

c c = -ig2C a - (k c /2)c c . (8) 

These equations describe both coherent evolution (terms proportional to £l(t),gi^) and population loss in the sys- 
tem (terms proportional to 7, k QjC ). The population loss in the system can be connected to direct radiative emission 
of the excited state |e) (at a rate 7|c e | 2 ), radiation leakage and absorption losses of mode a (K Q |c a | 2 ), and absorption 
and leakage out of mode c (k c |c c | 2 , of which K c ,ex\c c \ 2 is successfully out-coupled to a waveguide). In general the 
efficiency of extracting a single photon of frequency lo c out into the waveguide is thus 



J dt K c , ex \c c (t)\ 2 



J™dt Kc \c c (t)\ 2 +K a \c a (t)\2 +1 \c e (tW 



(9) 



For arbitrary f2(i), Eqs. ((8]) and ([9]) can be evaluated numerically. However, in certain limits one can find approximate 
solutions. In particular, when J7(i) and its rate of change are small compared to the natural oscillation and decay rates 
of the system, the state amplitudes c a , Cje will follow the instantaneous value of c s (t). Formally, we can adiabatically 
eliminate these states, setting Cj = for i — a,c, e. Then, one finds 

c s (t) = -^^c s (t), (10) 

7total 

while the other Ci(xc s (t), with the proportionality coefficients being functions of g±, g%, K a , k c , 0(f). The resulting 
substitution of the solutions of Ci(t) into Eq. ([9]) allows great simplification because the integrands now become 
time-independent, and after some simplification yields Eq. ([B]). Self-consistency of the adiabatic elimination solution 
requires that the the effective rate of population loss ~4il(i) 2 /7 to tai predicted from state \s) does not exceed the rate 
k c that a photon can leak out through the cavity mode c. 

In the effective wave- function approach used here, the population leakage out of mode c can also be explicitly related 
to the shape of the outgoing single-photon wavepacket. For instance, we can model the linear coupling of cavity mode 
c to photons propagating in a single direction in a waveguide with the following Hamiltonian (in a rotating frame), 

H w = / dkhv(k - uj c /v)a\a k - hg w dk {a\dke lkZc + h.c) . (11) 



Here k denotes the set of wavevectors of the continuum of waveguide modes, v is the velocity of waveguide fields, g w 
is the coupling strength between cavity and waveguide modes, and z c denotes the position along the waveguide where 
the cavity is coupled to it (for simplicity we set z c = from this point on). Since we are now explicitly accounting 
for the waveguide degrees of freedom, we add a term J dk Ck(t)\lk) to the effective wave- function of the system. The 
equations of motion of the total system are identical to Eq. © , except that 



c c = -ig%c a - (K Ctin /2)c c + ig w J dkc k , (12) 
c fe = -iv(Sk)c k + ig w c c , (13) 
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where 6k = k — u> c /v. Compared to Eq. |(5J), we have now included the coupling of mode c to the waveguide, and 
accordingly have replaced n c —>K c in in the equation for c c since the leakage into the waveguide should be accounted 
for by the new coupling terms. The equation for c k can be formally integrated; assuming that the waveguide initially 
is unoccupied, Cfc(O) = 0, one has 



c k (t)=ig w / dt / Cc (Oe~* cafc(t ~ t) . (14) 



Substituting this into the equation for c c and performing the Wigner-Weisskopf approximation 1411 , one recovers 
the expression for c c in Eq. |jHJ) by identifying K c . ex — 2ng^ u /v. The one-photon wave-function [4l| is given by 
ip w (z,t) — (vac\E w (z,t)\ip(0)) = (^/2nig w /v)0(z)c c (t — z/v), where &(z) is the step function. The wave-function 
shape is thus directly proportional to c c (t). Under adiabatic elimination, 

^(z, t) = ^^6(z) 8iglg2 2 M t z/v)c s {t - z/v), (15) 

and thus for a desired (and properly normalized) pulse shape vp w one needs only to solve Eqs. (p~5|) and (|10P to 
obtain the corresponding external field f2(f). It is straightforward to show that the normalization is given by 
J dz \ip w (z, i^oo)| 2 = F provided that c s (oo)— >0. This normalization reflects the probability that a single photon 
ends up in the waveguide. 



Nonlinear cavity design 

The nanobeam cavities are formed by a 4-period taper in the size and spacing of the holes in the uniform photonic 
"mirror" on both sides of the cavity center in order to introduce a localized potential for the TE and TM modes. The 
950-1425 nm cavity nanobeam has width w — 420 nm and depth d = 307.5 nm, and the hole spacing tapers from clq 
= 360 nm in the mirror to a c = 337 nm in the center. The holes were made elliptical to give an additional design 
parameter to separately optimize the TE and TM mode Q factors. The elliptical hole semi-axes are 84 nm and 108 
nm in the mirror section, and the hole size-to-spacing ratio is held constant through the taper section. This design 
yields cavity parameters of Q = 1.2 x 10 7 and V n = 0.77 for the TE mode, and Q = 7.3 x 10 4 and V n = 1.45 for 
the TM mode (V n is the mode volume normalized by (A/n) 3 ). The factor 7/70 = 0.10 (0.20) for the TE (TM) mode 
is determined by simulating the total power emitted by a non-resonant dipole source in the cavity center. We have 
accounted for the index dispersion of our candidate material, GaAs, for which n(1425 nm) = 3.38 and n(950 nm) = 

3.54 m. 

The nonlinear parameter gi is determined by calculating the volume integral of Eq. ([3]) using the exact mode fields, 

E a and E c , extracted from our 3D-FDTD calculation. Because the mode fields are oriented along the y and z-axes, 

respectively, as defined in Fig. [DJc), the classical field which drives the differency frequency generation, Ej,, must 

be polarized along x. This field has a frequency = uj a — lu c , which corresponds to a wavelength = 2.85 /^m. 

(2) (2) 
The relevant nonlinear susceptibility tensor elements are Xxyz(GaAs) = 2d\± = 550 pm/V and Xxyz(GaP) = 320 

pm/v mm. 

We assume the classical field is constant over the spatial extent of the cavity modes, which allows Eb tX to be taken 
in front of the integral for 52, giving 

92 = J dv X <£) z E: <y E c , z . (16) 

To justify this assumption, we simulated a Gaussian beam with \ = 2.85 /im that is focused by a lens with a modest 
numerical aperture (NA) of 0.5 onto a ridge waveguide, and found that the average field amplitude is approximately 
uniform over the linear extent of our cavity modes (approx. x = -1 fim to +1 /im). In the gi calculation, the 
magnitude of E^, x for a given beam power, Pb, is then determined from the relation = eocnr 2 E^ x / '4, where r is 
the focal spot radius. 
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